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Abstract. Electron-impact cross sections for the dissociation process of vibrationally 
excited He^ molecular ion, as a function of the incident electron energy are calculated 
for the dissociative transition —>■ A^I]+ by using the i?-matrix method in the 

adiabatic-nuclei approximation. The potential energy curves for the involved electronic 
states and target properties, also calculated with the i?-matrix method, were found to 
be in good agreement with the results reported in literature. 


1. Introduction 

He^ is a simple molecular ion that is present in many natural and technological 
systems. It plays a dominant role in electric discharges [Tj and stellar medium [2], 
as well as in the chemistry of the early universe [31 E] and thermonuclear fusion [5]. 
In helium containing plasmas, electron-impact dissociative recombination and direct 
dissociation are among the main reactive processes of the He^ ion leading to atomic 
and ionic species. Dissociative recombination is a low-energy process occurring through 
temporarily neutralization of He^ ion followed by dissociation into two He atoms. It 
has been widely studied due its importance in relatively low-temperature system (see 
Ref. [6] and references therein). On the opposite side, direct dissociation of He^ induced 
by electron-impact into He atom and He"*" ion requires fast electrons with a kinetic 
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energy equal or larger than the excitation energy of the process, which amounts, at the 
ground state equilibrium bond-length, to about 10 eV. Corresponding temperatures of 
this order of magnitude of energies are typical for the divertor plasmas of toroidal fusion 
devices. The plasma cooling in the divertor region favors the recombination processes 
which allow for the formation of molecular excited species, so that plasma chemistry, 
involving gas-phase and surface becomes rather complex. 

In the present work we undertake a theoretical study of cross sections for dissociative 
excitation induced by electron-impact which, to the best of our knowledge, has not been 
studied before theoretically or experimentally. The process can be represented as, 

e-^ Ret = 0-23) e+ Ret{A^^t) ^ e + He + He+, (1) 

where v represent the vibrational level of He^. The molecular states and 

are degenerate at inhnitely large internuclear distances R and the energy splitting of 
these two states at hnite i? is a result of the resonant one-electron exchange between 
the two ion cores He"*". 

The cross sections, calculated for all the 24 vibrational levels of the ground state 
potential well, were obtained by using the i?-matrix method |7j in connection with the 
adiabatic nuclei (AN) approximation [8|,[9]. The i?-matrix method, in fact, assumes an 
excitation of the target occurring at a hxed internuclear distance, at which the scattering 
T-matrix is calculated, so that the nuclear motion is suppressed and a purely electronic 
excitation cross section is produced. In order to take into account the vibrational 
dynamics of the process ([T]) in a more realistic way, we used in our calculations the 
well-known AN approximation, which is based on a vibrational average of the scattering 
T-matrix calculated for a sufficient number of internuclear distances. 

The paper is organized as follows: In Section [2] the i?-matrix method is outlined 
and the construction of target states within the hxed-nuclei approximation is discussed. 
In Section [3] the main equations of AN approximation are formulated and in Section 0] 
we present and discuss our results. Section [5l hnally, concludes the work. 

2. i?-matrix method 

The main facet of the i?-matrix method is the division of conhgurational space into an 
inner and an outer region. The inner region is dehned as the volume of a sphere of 
radius a centered at the center-of-mass of the target molecule. In this work we have 
set a = 12 a.u. This region is constructed so that the wave functions of all N except 
the scattering electron vanish at boundary of the sphere. In this region the exchange 
effects, i.e. the short range electron-electron correlations and polarization effects, are 
important. Implicitly this method assumes that the Pauli principle, which asserts that 
all inner region electrons are identical and any many electron wave function must be 
anti-symmetric to interchange of these electrons. In the outer region it is assumed that 
one electron can be considered to be distinct. This electron therefore moves in a local 
potential arising from its long-range interaction with the target. The exchange effects 
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are neglected in this region. The construction of matrix R provides the link between the 
inner and outer region. The i?-matrix method has some distinct advantages. The main 
advantage is that the inner region problem needs to be solved only once. The energy 
dependence is obtained entirely from the solution of the much simpler outer region case. 
This allows us to generate solutions at a large number of energies at a minimal extra 
computational effort. This generates cross sections on a hne energy grid that helps in 
analyzing the nature of resonances that may arise. 

In the inner region the total wave function is expanded in a conhguration-interaction 
(Cl) basis which takes the following form for each total orbital angular momentum, spin 
and parity combination: 

^ ^ (xi... xat) ^ji^N+i) aijk + ^ Xm(xi... xtv+i) bmk , (2) 

i,j m 

where A is an anti-symmetrization operator, xat are the spatial and spin coordinates of 
the iV-th electron, is the wave function of the Tth target state, are the continuum 
orbitals of the scattering electron, k represents a particular i?-matrix basis function. 
The variational coefficients a^jk and bmk are determined by matrix diagonalization. 
Furthermore the electrons, whose space-spin coordinates are represent by Xj, must obey 
the Pauli principle and are therefore anti-symmetrized by operator A. In practical 
implementations for generating the conhgurations which make up this term it is often 
necessary to impose a constraint on the coupling of the hrst N electrons to ensure 
that the target wave function does not get contaminated by states with the same 
conhguration but different space-spin symmetry. The second summation involves 
conhgurations which have no amplitude on the i?-matrix boundary and where all 
electrons are placed in orbitals associated with the target. Since they are conhned to a 
hnite volume of space they will be referred to as conhgurations. Such conhgurations 
are essential to relax the constraint of in orthogonalization between the continuum 
orbitals and those belonging to the target of the same symmetry. In more sophisticated 
models the conhgurations are also used to model the ehects of target polarization. 
The remaining information that is required to set-up the outer region problem concerns 
properties of the target. The target state energies relative to the ground state are needed 
as they give the energies of the asymptotic channels. The multipole moments associated 
with these target states are also required as they determine the outer region. We include 
only the dipole and the quadrupole transition moments in the present work. 

In our calculations we have 16 scattering symmetries. These are singlet and triplet 
spin states of the 8 irreducible representations of D 2 h point group. We have included 
up to g partial waves contribution to the cross sections, the ehect of the remaining 
partial waves is included via a Coulomb-Born closure approximation. The i?-matrix 
constructed at 12 a.u. is propagated up to 50 a.u. where the effect of multipole moments 
is negligible [lO] . The scattering solutions are then matched with the Coulomb boundary 
conditions to extract K matrices, T matrices and other observables. We have employed 
UKRMol code suite of programs m- 
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Figure 1: (a) Ground (X^S+) and excited (A^S+) i2-matrix electronic potential energies 
and (b) transition dipole moments compared with the calculations taken from Refs. [12] and 
m respectively. 

2.1. Target states 

At the i?-matrix boundary between the inner and outer regions, the amplitude of the 
target orbitals must be negligible. This restricts the use of extended diffuse functions in 
a chosen basis set. Also a single basis set is mandatory to represent all the target states 
included in the calculation. 

The basis set employed in this work is the cc-pVTZ Gaussian basis set for He^ 
molecule. This set includes polarization functions. The electronic conhguration for the 
ground state of He^ in its natural symmetry is (lug)^ (lo'u)- The molecule is 

treated in a reduced D 2 h symmetry in which there are eight symmetries: Ag, i?ig, 
Biu, B 2 g, Bsg, i? 2 u, B^u- The calculations were performed at the different equilibrium 
bond-lengths required in the calculation of AN cross sections. In particular for the 
equilibrium distance of 2.042 a.n. the self-consistent held (SCF) energy for the X state 
is —4.92150 a.n. In the close conpling expansion of the trial wave function of the 
scattering system, we included the X state and the hrst excited state A^S+, {ICgY 
(2(Tg). Cl wave functions are used to represent each target state. In our Cl model, 
we have the occupied orbitals which are augmented by the virtual (vacant) molecular 
orbitals. We include virtual orbitals up to Ta^, 362«, ^bsu, 362g, ^b^g, la^. 

The three molecular orbitals were free to move into the entire active space. The vertical 
excitation energy of the A state at equilibrium distance is 9.954 eV with respect to the 
X state. We obtained for the dipole transition moment of X-A transition a value of 
0.96856 a.u. 

Figures [I](a) and {b) show a comparison of our i?-matrix potential curves and 
transition dipole moments, with those of Refs. [12] and [13] respectively. The agreement 
is quite good in both cases. The ground and excited states of He^ support respectively 
24 and 2 vibrational levels. The corresponding eigenvalues, calculated by using the 
curves of Ref. [12], are shown in Table [H The shallow minimum of A state appears at 
R ~ 8.7 a.u. an internuclear distance far away from the Franck-Condon (FC) regions 
of X so it cannot affect the dissociation process in any way. 
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H4 


V 

Present results 

Ref. [H] 

Present results 

Ref. 03 

0 

0.1042 

0.1042 

2.4733 

2.4730 

1 

0.3069 

0.3061 

2.4743 

2.4740 

2 

0.4997 

0.4992 



3 

0.6839 

0.6837 



4 

0.8597 

0.8594 



5 

1.0265 

1.0264 



6 

1.1846 

1.1846 



7 

1.3346 

1.3340 



8 

1.4761 

1.4745 



9 

1.6094 

1.6062 



10 

1.7347 

1.7289 



11 

1.8521 

1.8426 



12 

1.9607 

1.9472 



13 

2.0597 

2.0426 



14 

2.1485 

2.1288 



15 

2.2266 

2.2056 



16 

2.2938 

2.2729 



17 

2.3501 

2.3306 



18 

2.3952 

2.3785 



19 

2.4292 

2.4165 



20 

2.4523 

2.4446 



21 

2.4651 

2.4627 



22 

2.4700 

2.4715 



23 

2.4740 

2.4739 




Table 1: Vibrational eigenvalues (eV) for the and electronic states compared 

with those from Ref. m and Ref. [12] respectively. All the values are referred to the bottom 
of the ground state potential curve. 


3. Adiabatic nuclei approximation 

The AN approximation is well-known in scattering theory and the interested reader can 
refer to the rich literature on the subject [HI El for details. We will only recall here the 
main aspects of its formulation for diatomic molecules and state the relevant equations 
useful for our discussion. 

The AN in its standard treatment, is based on the decoupling of the nuclear and 
electronic motion, so that the total scattering wave function can be approximated 
by the product of the electronic wave function in the hxed-nuclei approximation, i.e. 
parametrically dependent on the internuclear distance R, and the nuclear ro-vibrational 
wave function. This factorization is based on the assumption that the speed of the 
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target and incident electrons is much higher than that of the nuclei, so that an electronic 
transitions take place in a very short time, during which the nuclei remain still at a given 
internuclear distance. In this paper we are interested on vibrationally excited target and 
then we limit the calculation to rotationally averaged cross sections only [8]. Moreover, 
as the final electronic state A is essentially repulsive, so that its vibrational spectrum 
form a continuum and the excitation to one of these vibrational levels is followed by 
dissociation, the final total dissociation cross section implies an integration over the 
continuum vibrational spectrum. Namely: 

dCr 


cr„ 


^X,A/ 


[e = 


^th 


der_ 


TT 


9s dCc 
'.Ml' 9 tth 




T,y(R; k) 


■ ( 3 ) 


where e is the incident electron energy and k = (2 rUe its momentum. k) 

is the /^-depending T-matrix expressed in terms of the I, V partial wave quantum 
numbers, A is the group symmetry index and S is the spin state quantum number 
of the global electron-molecule system, gs is a spin multiplicity factor, given as 
(25' -I- l)/2(2S'j -|- 1), where 2(25'* -|- 1) is the initial total spin. It assumes the values 1/4 
and 3/4. I and V run over the values from 1 to 3 and 1 to 9 respectively, 1 < A < 16 for 
the D 2 h symmetry group. Cc is the continuum energy and Cth the continuum threshold. 
Finally, E is the total energy dehned by i? = e -|- e„, where Cy is the initial vibrational 
energy eigenvalue. The T-matrix has been calculated with the i?-matrix method for 
1000 values of the internuclear distance R in the interval [1.0, 6.0] a.u. for 0 < n < 15 
and for 1800 i?-values in the interval [1.0,10.0] a.u. for v > 15. This is done because 
the continuum wave functions extend to large internuclear distances. 

Eq. ([3]) can be cast in an approximate form. Assuming, in fact. 




dCr 


'^th 


dCr 


( 4 ) 


using the closure property J dcc |c)(c| = 1 and applying then the FC approximation at 
a particular value of the internuclear distance R = R, it can be written as: 

2 


(T,. 




[e = 


TT 


ds 




Tff(R\ k) 


1 , 1 ' 


( 5 ) 


This is the approximate form of the total cross section generally used in i?—matrix 
calculations, where, usually, R is taken at the equilibrium internuclear distance. It will 
be useful in the discussion of our results in the following sections. 


4. Results 

Figure [2][a) shows the cross sections for n = 0, as a function of the incident electron 
energy, calculated in the AN approximation by Eq. (|3]) (full-blue line). The cross section 
curve exhibits a smooth behavior in all the considered range of energies, reaches a 
maximum at 14-15 eV and then decreases for higher energies. A structure, characterized 
by a small maximum positioned at about 7 eV, appears in the energy range of 6.5-8 eV. 
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Electron energy (eV) 


Internuclear distance (a.u.) 


Figure 2: (o) Cross section as a function of the incident electron energy for process © 
starting from r = 0, calculated in the AN (full-blue line) and in the fixed-nuclei (dot-red line) 
approximations. (6) Potential energy curves and u = 0 wave function (arbitrary units) as a 
function of the internuclear distance. The two vertical dashed lines, starting from the classical 
turning points of the zero-level, enclose the FC region. The zero-energy is placed at the bottom 
of the ground state well. 


The general trend of the cross section can be explained by inspecting the behavior of 
the vibrational wave function for the level u = 0 in connection with the electronic state 
potential curves as depicted in Fig. |2](6). The right panel of the hgure shows that the 
threshold of the process is at 2.369 eV, counted from the u = 0 level, but the apparent 
threshold is observed at ~ 6.5 eV (see Fig. [2](a)), corresponding to the transition energy 
a.t R = 2.39 a.u., close to the outer limit of the FC region. 

The local maximum at 7 eV is generated by the formation of a resonance occurring 
in the low energy region. This can be better appreciated in Fig. [3l where the cross 
sections calculated by Eq. ([5]) for a number of bond-lengths, are shown. The hxed- 
nuclei cross sections calculated at i? = 2.39 a.u. and R = 2.33 a.u. exhibits a sharp 
peak close to 7 eV, while in the other curves, calculated for decreasing internuclear 
distances, the maximum gets reduced and hnally disappears for R > 2.13 a.u. This is 
consistent with fact that the threshold of the process, correspondingly, increases and 
at large collision energies the formation of a resonant state is prevented by the small 
probability of electron capture. The resonance maximum observed at the two largest 
bond-lengths in Fig. [3l induces the formation of the 7-eV peak in the AN cross sections, 
whose magnitude, however, is greatly reduced with respect to the hxed-nuclei cross 
section. This is due to the fact that the two above resonance bond-length values fall in 
the right edge of the FC region, where the n = 0 vibrational wave function is going to 
vanish (see Fig. [2](b)). 

In Figure H] the AN cross sections are shown for different vibrational levels. For 
these cases, the peak maximum increases up to n = 2, then reduces and its position 
gets slightly shifted toward lower energies. This behavior is probably due again to the 
modulation of the vibrational wave function which, for n > 0, is peaked at the right edge 
of the FC region and this could enhance the resonance peak arising in the hxed-nuclei 
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Figure 3: i2-matrix cross section calculated at the bond-lengths shown in the plot. Reg 
denotes the equilibrium internuclear distance. 



Figure 4: AN cross sections for processes starting from different He^ vibrational levels v. 

cross sections. For higher levels, instead, the FC region extend far from the resonance 
bond-lengths so that the wave function maximum no longer overlap with the resonance 
peak and the AN cross section maximum becomes less pronounced. 

At high energies, the AN cross section for v = 0, in Fig. |2]^a), strictly follows 
the hxed-nuclei cross section calculated at the equilibrium internuclear distance, 
Reg = 2.042 a.u. (red-dotted curve), showing that for high energies the hxed-nuclei 
approximation furnishes quite accurate results. For higher vibrational levels. Fig. 0] 
shows a decreasing trend of the cross section for large energies. This is a counter intuitive 
behavior if one considers that the threshold energies gets smaller with the increasing 
of V. Probably, again the vibrational wave function, whose amplitude decreases for 
high levels, plays a decisive role in affecting the cross sections. Finally, the oscillations 
observed in the cross section curves for n > 0 can be attributed to oscillating behavior 
of the FC densities also observed in other circumstances na. 

5. Conclusions 


We have calculated the adiabatic nuclei approximation cross sections for electron-impact 
dissociation of vibrationally excited He^ cation. For the case n = 0, the cross section has 
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a smooth energy behavior, presenting a structure only at the low energy generated by 
the formation of a resonant state. The intensity of the resonance peak is reduced by the 
small values of the decaying part of the vibrational wave function in the external side 
of the FC region, where the resonance arises. For n > 0 the structure is enhanced, 
with respect to the previous case, by the corresponding vibrational wave functions 
which reach their maximum at the outer edge of the FC region. At larger energies, the 
general decreasing trend of the cross sections is likely determined, once again, by the 
corresponding behavior of the initial state vibrational wave function amplitude, while 
the oscillating structures of the cross section curves can be attributed to the influence 
of the continuum FC density. 
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